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We present a formalism for obtaining the statistical properties of functionals and inverse functionals of the 
paths of a particle diffusing in a one-dimensional quenched random potential. We demonstrate the implemen- 
tation of the formalism in two specific examples: (1) where the functional corresponds to the local time spent 
by the particle around the origin and (2) where the functional corresponds to the occupation time spent by the 
particle on the positive side of the origin, within an observation time window of size t. We compute the disorder 
average distributions of the local time, the inverse local time, the occupation time and the inverse occupation 
time, and show that in many cases disorder modifies the behavior drastically. 
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I. INTRODUCTION 

The statistical properties of functionals of a one dimen- 
sional Brownian motion have been extensively studied and 
have found numerous applications in diverse fields ranging 
from probability theory yj, i3j], finance |5l meso- 
copic physics 01, computer science j^l, and in understanding 
weather records f9'l. The position x{t) of a one dimensional 
Brownian motion evolves with time t via the Langevin equa- 
tion 



starting from a;(0) = xq, where ?7(t) is a thermal Gaus- 
sian white noise with mean (77 (r)) = and a correlator 
(77(t)?7(t')) = 5(t — t'). a functional T is simply the in- 
tegral up to time t 

T= f V{x{T))dT, (2) 
' Jo 

where V{x) is a prescribed non-negative function whose 
choice depends on the specific observable of interest. For a 
fixed initial position of the Brownian motion and a fixed 
observation time t, the value of T varies from one history or 
realization of the Brownian path {x{t)} to another (see Fig.Q 
and a natural question is: what is the probability density func- 
tion (pdf) P(r|t, xo)7 
[ Following the path integral methods devised by Feyn- 
man fioll . Kac showed (J 01 that the calculation of the 
pdf P{T\t, xo) can essentially be reduced to a quantum me- 
chanics problem, namely solving a single particle Shrodinger 
equation in an external potential V{x). This formalism is 
known in the literature as the celebrated Feynman-Kac for- 
mula. Subsequently, this method has been widely used to 
calculate the pdf of T with different choices of V{x) as de- 
manded by specific applications. This has been reviewed re- 
cently in Ref. fs'l. In particular, two most popular applications 
coiTespond respectively to the choices, V{x) = 6{x — a) 
and V{x) — 9{x), where 5{x) is the Dirac's delta func- 
tion and 9{x) is the Heaviside step function. In the former 




FIG. 1: Schematic plots of T defined by Eq. J2j as a func- 
tion of t, corresponding to four different realizations of the paths 
[{a;(r)}, for < r < denoted by Ri,R2,Rs and R4 respec- 
tively. For fixed t (ti or t2 or ts or t4, shown by vertical dashed 
lines), T takes different value for different realizations. On the other 
hand for a fixed T (horizontal dashed line) corresponding t is differ- 
ent for different realizations: ti for Ri, t2 for R2, ta for R3 and t4 
for R4. 



case, the corresponding functional T{a) = S{x{t) — a) dr 
has the following physical meaning: T(a) da is just the time 
spent by the particle in the vicinity of the point a in space, 
i.e., in the region [a, a + da], out of the total observation 
time t. Note that, by definition, jT{a)da — t. The func- 
tional T{a) is known as the 'local time' (density) in the lit- 
erature. In the second case V{x) — 9{x), the functional 
T = 9{x{t)) dr measures the time spent by the particle 
on the positive side of the origin out of the total time t and is 
known as the 'occupation' time. The probability distribution 
of the occupation time was originally computed by Levy lil ill . 
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/(f P{T'\t, 0) dT' = f arcsin (^57^)' ™d is known as the 
arcsine law of Levy. Since then, the local and the occupa- 
tion times for pure diffusion, have been studied extensively 
by mathematicians (Illll H d d d. Recently, the 
study of the occupation time has seen a revival in physics 
literature and has been used in understanding the dynamics 
out of equilibrium in coarsening systems |18, 19], ergodicity 
properties in anomalously diffusive processes L20n21il . in re- 
newal processes i22il . in models related to s pin glasses i23ll . 
in simple models of blinking quantum dots |24], and also in 
the context of persistence L25.. .26.1 . Local and Occupation 
times have been also studied in the context of stochastic er- 
godicity breaking |27], first passage time |28], diffusion con- 
trolled reactions activated by catalytic sites |29j and diffusion 
on graphs HsollsTll . In polymer science, a long flexible poly- 
mer of length t is often modeled by a Brownian path up to 
time t. In this context, the local time at a position j^is propor- 
tional to the concentration of monomers at r in a polymer of 
length t. 

A natural and important question is how to generalize the 
Feynman-Kac formalism to calculate the statistical properties 
of the functionals of the type in Eq. when x{t) is not just 
a pure diffusion process, but it represents the position of a 
particle in an external random medium. While various prop- 
erties of diffusion in random media have been widely studied 
in the past ll32ll33ll34ll35ll . the study of statistical properties 
of functionals in random media is yet to receive its much de- 
served attention. In this paper we undertake this task. More 
precisely, we are interested in calculating the pdf P{T\t, xo) 
of a functional T as in Eq. (|2} where x{t) now evolves via the 
Langevin equation 

^^F{x{r))+Tj{r) (3) 

where 7/(t) represents the thermal noise as in Eq. ([0 and 
F{x) = —dU /dx represents the external force, the deriva- 
tive of the potential U (x), felt by the particle. Most generally, 
the external potential consists of two parts, U{x) = Ud(x) + 
Ur{x), a deterministic part Ud{x) and a random part Ur{x). 
The random part of the potential Ur{x) is 'quenched' in the 
sense that it does not change during the time evolution of the 
particle, but fluctuates from one sample to another according 
to some prescribed probability distribution. Consequently, the 
pdf P{T\t, xo) will also fluctuate from one sample of the ran- 
dom potential to another and the goal is to compute the disor- 
der averaged pdf P{T\t,xo) where the — denote the average 
over the distribution of the random potential. A popular model 
for the random potential is the celebrated Sinai model L36.I . 
where various disorder averaged physical quantities can be 
computed analytically fslEilliSllliSEllEl, and 
yet the results exhibit rich and nontrivial behaviors and also 
capture many of the qualitative behaviors of more complex 
realistic disordered systems. The Sinai model assumes that 
Ur{x) = -JaB{x) where B{x) represents a Brownian motion 
in space, i.e.. 



where ^{x) is a Gaussian white noise with mean {£,{x)) = 
and a correlator {^{x)£^{x')) = 5{x — x'). The constant a 
represents the strength of the random potential. 

In this paper, we first present a generalization of the 
Feynman-Kac formalism to calculate the pdf P(r|t,a;o) in 
presence of an arbitrary external potential U{x). To obtain 
explicit results using this formalism, we next assume that the 
random part of the potential is as in the Sinai model, i.e., our 
external potential is of the form U{x) — Ud{x) + y/aB{x), 
where B{x) is a Brownian motion in space and Ud{x) is 
the non-random deterministic part of the potential. It turns 
out that the asymptotic behavior of the disorder averaged 
pdf P{T\t,XQ), quite generically, has three different quali- 
tative behavior depending on the curvature of the determin- 
istic potential Ud{x), i.e., whether Ud{x) has a convex (con- 
cave upward) shape representing a stable potential (i.e., attrac- 
tive force towards the origin), a concave (concave downward) 
shape representing unstable potential (a repulsive force away 
from the origin) or just flat indicating the absence of any ex- 
ternal potential. To facilitate explicit calculation, we model 
the deterministic potential simply by, Ud{x) — ~n\x\, so that 
/i < represents a stable potential, /i > represents an un- 
stable potential and = represents a flat potential. This 
specific choice facilitates explicit calculation, but the results 
are qualitatively similar if one chooses another form of this 
potential. Thus, in our model, we will consider the external 
potential as 

U{x) = -^fi\x\ + V^B{x) (5) 

where B{x) = ^{x') dx' is the trajectory of a Brownian 
motion in space (see Fig. |2j. Note that the case — cor- 
responds to the pure Sinai model. Figure (|2jl shows typical 
realization of potentials for /i = 0, /i > and /i < 0. The 
corresponding force in Eq. (|3} is simply given by 

F(x) = fisign{x) + y/a£_{x). (6) 

We will demonstrate how to calculate exphcitly, using our 
generalized Feynman-Kac formalism, the disorder averaged 
pdf P(r|t, xq) when the external potential is of the form given 
by Eq. (|5jl. Despite the simplicity of the choice of the exter- 
nal potential, a variety of rich and interesting behavior can be 
obtained by tuning the parameter ji/cr, as shown in this pa- 
per We will present detailed results for the two functionals, 
namely for the local time and the occupation time correspond- 
ing to the choices V{x) ~ 5{x) and V{x) — 9{x) respectively 
in Eq. (|2jl. Also, to keep the discussion simple, we will present 
our final results for xo = corresponding to the particle start- 
ing at the origin. However, our method is not limited only to 
this specific choice. Some of these results were briefly an- 
nounced in a previous letter . 

In addition, in this paper we also introduce the notion of 
'inverse functional', which is defined as follows. If V{x) in 
Eq. (|2jl is non-negative, then for each path {x{t)}, T is a 
non-decreasing function of t, which we formally denote by 
T = g(t|{a;(r)}, xq). Therefore for a given realization of 
path {a;(r)} and given T there is a unique value of t (see 
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FIG. 2: A classical particle (represented by •) diffusing in a typical realization of the potential U{x) = —fi\x\ + ^/aB{x), where B{x) 
represents the trajectory of a Brownian motion in space with B{0) — 0. The three figures are for /i = 0, /i > and /i < respectively. The 
dash lines show the potential for a = 0. 



Fig. n), which we formally write as the inverse of the func- 
tional g i49il 



t = g-'(T\{x{r)},xo). 



(7) 



This inverse time t physically means the observation time that 
is required for any given path {x{t)} in order to produce a 
prescribed value of T. Of course, for the same value T, for a 
different path {x{t)}, the value of t will be different. Thus, t 
is a random variable for a fixed T, which takes different values 
for different realizations of paths and we would like to com- 
pute its pdf, which we denote by I{t\T,xo) and by definition 
I{t\T, 0) dt = 1. Clearly, this pdf will also differ from 
sample to sample of the external potential in Eq. (|5Jl and our 
goal is to obtain the disorder averaged distribution I{t\T, xq). 
In this paper, we present detailed results for I{t\T,Q) again 
for the two choices of V{x) = 5{x) and V{x) = 9{x) cor- 
responding to the local time and the occupation time respec- 
tively. The inverse local and occupation times are useful for 
experimentalists as they provide an estimate on the required 
measurement time. For example, in the context of polymers, 
the inverse local time is the typical length of a polymer re- 
quired to obtain a certain monomer concentration. 

The rest of the paper is organized as follows. In SecHH we 
present our general approach for computing the pdf P(T|t, x) 
of the functional T defined by Eq. (|3 for a given t, and the 
pdf /(t|T, x) of the inverse functional defined by Eq. for 
a given T, for a given sample of the random potential, for 
arbitrary starting position of the particle a;(0) = x and for 
arbitrary but non-negative V{x). After this section we con- 
sider the specific examples of local time and occupation time 
by setting V{x) — 5{x) and V{x) = 6{x) respectively. We 
will use different notations for the pdfs in the two examples to 
avoid any misunderstanding. In the first example, where T is 
the local time, we denote the pdf of the local time P{T\t, 0) 
for a given t by P\oc{T\t), and the pdf of the inverse local 
time I{t\T, 0) for a given T by /ioc(^|r). In the second ex- 
ample, where T is the occupation time we denote the pdf of 
the occupation time P{T\t, 0) for a given t by Pocc{T\t) and 
the pdf of the inverse occupation time /(ijT, 0) for a given 
T by /occ(i|r). While our final goal is to obtain the disor- 
der averaged distributions P\oc{T\t), I\oc{t\T), Pocc(T|t) and 
-focc(i|r), it is however, instructive to study the pure case first. 



before tackling the problem with disorder which is obviously 
harder. In the same spirit, we have presented the detailed 
discussions on the local time, inverse local time, occupation 
time and inverse occupation time for the pure case (cr = 0) 
in Sees. |lll] |V] IVIII and IIXI respectivelv. before computing 
their disorder average in Sees . II VI I VII I Vllll andlXlrespectivel v. 
Sec. lXII contains some concluding remarks. Some of the de- 
tails are relegated to the appendixes. The results are summa- 
rized in TableslHllll andCni 



II. GENERAL APPROACH 

In this section we will show how to compute the pdfs 
P{T\t,x) and I{t\T,x) for arbitrary non-negative V{x) and 
arbitrary starting position x(0) — x, for each realization 
of random force F{x), by using a backward Fokker-Planck 
equation approach. In the following discussion we will de- 
note the functional defined in Eq. (|2j by g(t\{x{T)},xo), 
and use T as the value of the functional for a given path 
[{x{t)}, forO < t < t]. 

Since V{x) is considered to be non-negative, T defined by 
Eq. (|2} has only positive support. Therefore, a natural step is 
to introduce the Laplace transform of the pdf P{T\t, x) with 
respect to T 



P{T\t,x)e-P'^ dT 
-pg{t\{x(T)} ,x) 



exp 



x{id)=x 



V[x{t')] dt' 



(8) 



x{id)=x 



where 



lx{0)-- 



. denotes the average over all paths that start at 
the position a;(0) = x and propagate up to time t. Our aim 
is to derive a backward Fokker-Planck equation for Qp{x, t) 
with respect to the initial position x{0) — x. 

We consider a particle starting from the initial position x, 
evolves via Eq. Q up to time t + At. Now we split the 
time interval [0, t + At] into two parts: an infinitesimal in- 
terval [0, At], over which the particle experiences an infinites- 
imal displacement Ax from its initial position x and the re- 
maining interval [At, t + At] in which the particle evolves 
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from a starting position x + Ax. Since x(0) = x, one gets 

/J^V[x(t')]rf^' = V{x)At + 0[{Atf]. Therefore using 
Eq. (|8}, and splitting the integral over t' into the above two 
time intervals we obtain 



Qp{x,t + At) ^ (cxp l-p J 



V[x{t')]dt' 



x{0)=x 



= e-^'^(^)'^*(gp(x + Ax,t)\ , (9) 

where ( )ax denotes the average with respect to all possible 
displacements Ax. 

Now, in the limit At 0, integrating Eq. (|3} one gets, 

pAt 

Ax^F{x)At+ r]{T) dr + 0[{Atf]. (10) 

^0 

Hence, using the zero mean and the uncorrected properties of 
the noise we get 



limiM 

At^O At 



F(x) and lim 
At^a 



At 



1, 



(11) 



Therefore, by Taylor expanding Eq. (|9} for small Ax, and tak- 
ing the limit At 0, one arrives at the 'backward' Fokker- 
Planck equation. 



dQp 
dt 



2 9x2 



Fix) 



dQp 
dx 



pV{x)Qp, (12) 



with the initial condition Qp{x,0) = 1, which is easily 
checked by Eq. (|8}. The advantage of the above equation 
over the usual Feynman-Kac formalism 1 1, 2] is that, in the 
later case one has a 'forward' Fokker-Planck equation (spa- 
tial derivative with respect to the final position) where after 
obtaining the solution of the differential equation, one has to 
again perform an additional step of integration over the final 
position. In contrast, Eq. (I12> involves the spatial derivatives 
with respect to the initial position of the particle, and hence 
no additional step of integration over the final position is re- 
quired. 

The standard practice of attacking the partial differential 
equations of above type is by using the method of Laplace 
transform. We define the Laplace transform of Qp{x, t) with 
respect to t 



m(x) 



Qp{x,t)e-°'^ dt 



/ dte-"' / dTe-P^P{T\t,x), (13) 
Jo Jq 



where for notational convenience, we have suppressed the a 
and p dependence of u{x). Now by taking Laplace transform 
of Eq. il2i with respect to t we obtain the ordinary differential 
equation 
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u"{x) + F{x)u'{x) - [a+pV{x)]u{x) = -1, (14) 



where u'{x) — du/dx. The appropriate boundary conditions 
m(x — > ±cx)) are to be derived from the observation that if 



the particle starts at x ±00 it will never cross the origin 
in finite time. Note that Eq. J14> is valid for each sample of 
the quenched random force F{x). Thus in principle, from the 
solution u{x) one obtains P{T\t,x) by inverting the double 
Laplace transform in Eq. ( I13> for each sample of quenched 
random potential, and then takes the average over the disorder 
Our next goal in this section, is to show how to compute the 
pdf I{t\T, x) for a given sample of the quenched random force 
F(x). It turns out that I{t\T, x) is related to the pdf P{T\t, x) 
in their Laplace space as shown below. By definition we have. 



I{t\T,x)^{6{t-g-\T\{x{T)},x)))^ 



(15) 



However, it is elementary that for each realization of path 

{x(r)} 



dT 



5(t-g-'{T\{x{T)},x))=s{T~g{t\{x{T)},x)) 

(16) 

where \dT/dt\ is the usual Jacobian of the transformation, 
which is simply dT/dt as both T and t have only positive 
support. It immediately follows from the above two equations 
that 

Iit\T,x)^(^6(T-g{t\{x{r)},x))^^ . (17) 

Therefore, Laplace transform of I{t\T,x) with respect to T 
reads 

dTe-P^I{t\T,x) = (e-Ps{t\{'Ar)}.x) 



pdl' 



}p{x, t) , 



(18) 



where Qp{x,t) is given by Eq. (|8}. Now taking a further 
Laplace transform in Eq. ( I18> with respect to t, it is straight- 
forward to obtain 



/>00 /"OO 

/ die""*/ dTe~P^I{t\T,x) ^ 
Jo Jo 



1 — au{x) 



(19) 



Thus, we have established via Eqs. (I18> and ( I19> , the relation- 
ships between the Laplace transforms of the pdf of the func- 
tional T defined by Eq. (|2j and the pdf of the inverse func- 
tional defined by Eq. 0. Hence, again in principle, from the 
solution m(x) of the ordinary differential equation (I14> . one 
obtains I{t\T, x) by inverting the double Laplace transform in 
Eq. jl9> for each sample of quenched random potential, and 
then takes the average over the disorder Note that putting 
a = in Eq. ( I19> and inverting the Laplace transform with 
respect to p immediately verifies the normaUzation condition 
I{t\T,x) dt = 1. 
In the rest of the paper, we will demonstrate how to imple- 
ment this formalism for the particular examples of the local 
time corresponding to the choice V{x) = 6{x) and the occu- 
pation time corresponding to the choice V{x) = 9{x). Since 
in these examples we consider the starting position of the par- 
ticle to be the origin, we need to only find the solution u{0) of 
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the differential equation (I14> . In each example, we will con- 
sider the pure cases [a = 0) first, which help us anticipate the 
general features of the results in the disordered case (a > 0) 
studied later 



III. LOCAL TIME WITHOUT DISORDER (a = 0) 

In this case V{x) — S{x), corresponds to T in Eq. (|3 being 
the local time in the vicinity of the origin and P{T\t,0) in 
Eq. O being Pioc{T\t) — the pdf of the local time T for 
a given observation time window of size t and the starting 
position of the particle a;(0) = 0. For our purpose we only 
need the solution u{0) of the differential equation (I14> . which 
corresponds to the starting position of the particle being the 
origin. However, to obtain u{0) we have to solve Eq. il4\ in 
the entire region of x with the boundary conditions u{x 
zLoo) which are derived from the following observation. If 
the initial position x — > ±cxd, the particle can not reach the 
origin in finite time, which means that the local time T = 
0. Therefore, by substituting P{T\t,x ±oo) S{T) in 
Eq. ( I13> one obtains the boundary conditions 



u{x ±cxd) 



1 



Pio,(T|t), as 

poo poo 

Jo Jq 

X{a) 



a[p + A (a)] 
where X{a) is simply given by 



(25) 



z_(0)-z+(0) . y'±{x) 

= with z±{x) = — — . (26) 

2 y±{x) 

Note that putting p = in Eq. ( I25l l and then inverting the 
Laplace transform with respect to a readily verifies the nor- 
malization condition 



/>oo 

/ Pioc(T|t)dT = l. 

jQ 



(27) 



(20) 



Since A(a) is independent of p, inverting the Laplace trans- 
form in Eq. i25\ with respect to p yields 



(28) 



G(a) = / dte-"*P,oc(T|0 
Jo 

= ^ exp [-A(a)T] , 



We have to obtain the solutions u{x) = u+{x) for a; > and 
u{x) = u- (x) for X < Ohy solving Eq. J14t separately in the 
respective two regions. 
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u'^{x) + F(x)u'j.{x) — au±{x) = — 1, 



(21) 



with the boundary conditions u+{x oo) — 1/a and 
u-{x — oo) = 1/a, and then matching the two solutions 
u+ (x) and u- (x) at a; = 0. The matching conditions are 

u+(0) = M_(0) = u(0), u+(0) -u'_(0) = 2pu(0). (22) 

The first condition follows from the fact that the solution must 
be continuous at cc = and the second one is derived by inte- 
grating Eq. (I14> across a; = 0. 

By making a constant shift ii± (x) = l/a + A±y±{x),fiom 
Eq. i21\ one finds that y± (x) satisfy the p-independent homo- 
geneous equation 



^y'lix) + F{x)y'^{x) 



ay± (x) = 



(23) 



with the boundary conditions y+{x oo) — s- and y-{x 
— oo) — > 0. The constants A± are determined by the matching 
conditions given in Eq. (I22> . which can be rewritten as 



which is valid for any arbitrary force F{x). In the following 
subsections we will consider qualitatively different types of 
deterministic potentials to derive more explicit results. 



A. Flat potential 

We first consider the simple Brownian motion without any 
external force, F{x) = 0. In this case the solutions of Eq. i23i 
are obtained as 



y±ix) = 2/±(0)exp 



(29) 



Using the solutions in Eq. i26\ one gets A (a) = \/2a and 
hence the Laplace transform G{a) in Eq. ( I28> becomes 



G(a) = / dte-"*Pioc(T|t) 
Jo 



V2 



/2aT 



(30) 



Now inverting the Laplace transform with respect to a, one 
finds that the distribution of the local time is Gaussian for all 
T and t. 



PiociT\t) - ^exp 



nt 



' 2t 



(31) 



A+y+{0) = A_y_{0)=u{0)--, (24a) 

a 

A+j/V(0) - A^y'_iO) = 2pu{0). (24b) 

Eliminating the constants A± from Eq. (I24> . we obtain the 
Laplace transform w(0), defined by Eq. il3\ with P{T\t, 0) = 



B. Unstable potential 

Now we consider the case of a Brownian particle moving in 
an unstable potential U{x) such that U{x ±oo) — oo. 
The corresponding repulsive force F{x) drives the particle 
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eventually either to +00 or to —00. The pdf of the local time 
Pioc{T\t) in the case of an unstable potential tends to a steady 
distribution Pioc(T') as t ^ 00, which can be computed ex- 
plicitly. To see this consider the function G{a) in Eq. (I28> . By 
making a change of variable r = at, it follows from Eq. ( I28> . 



G(a) 



a JO 



(32) 



Assuming Pioc{T\t — > 00) = Pioc{T), we find form the above 
equation that G{a) P\ociT)/a as a ^ 0. Comparing this 
behavior with Eq. (l28t gives 

i^ioc(r) - A(0)exphA(0)T], (33) 

provided A(0) is a finite positive number Thus generic ally 
for all repulsive force F{x), the local time distribution has a 
universal Poisson distribution in the limit t ^ 00. The only 
dependence on the precise form of the force F{x) is through 
the rate constant A(0). 

The rate constant A(0) can be expressed in terms of the 
force F{x) in a more explicit manner. Putting a = in 
Eq. (I23> and solving the resulting equation with the bound- 
ary conditions y+{x 00) and y„ {x — * —00) we 
get 



y+ix) = y+{0) 



x > 0, 



(34) 



ih'^(y)dy 

y^{x)=y^{O y-^ \ x<0, (35) 

LaoV yy) dy 

where ■0(2/) = exp[— F{x) dx\. Substituting these results 
in Eq. M6\ gives the rate constant as 



A(0) 



t^i^Hy)dy !^r{y)dy^ 



(36) 



Let us now consider a simple example where the potential 
U{x) = — mI^^I with /i > 0, corresponding to the repulsive 
force F{x) ~ /isign(a;) from the origin. In this case ?A(y) = 
exp[— and hence from Eq. ( l36l we get A(0) = 2^. 



C. Stable potential 

We now turn our attention to the complementary situation 
when the potential U (x) is stable, i.e., U{x ±00) 00. In 
this case the force F{x) is attractive towards the origin so that 
the system eventually reaches a well defined stationary state. 
The stationary probability distribution p{x) for the position of 
the particle is given by the Gibbs measure 



p{x) 



-2U{x) 



(37) 



where U{x) — — F{x') dx' and Z is the partition function. 



Z 



(38) 



In this case the Laplace transform G{a) of the pdf of the lo- 
cal time Pioc{T\t) is still given by Eq. ( I28t . However, unlike 
the unstable potential in the previous section, the distribution 
Pioc{T\t) does not approach a steady state ast^oo. Instead 
it has a rather different asymptotic behavior. 

To deduce this asymptotic behavior, let us first consider the 
average local time (T) = {S[x{t')]) dt' . For large t' , the av- 
erage {5[x{t')]) approaches its stationary value {&[x{t')]) 
p(0), where p(0) — 1/Z from Eq. Hence ast —t 00 the 
ratio T /t approaches the limit 



t 



1 



(39) 



where Z is given by Eq. ( I38l l. Thus for large t, the average 
local time scales linearly with time t, which indicates that the 
natural scaling limit in this case is when t — > oo,T ^ 00 
but keeping the ratio r ~ T/t fixed. We will see that in this 
scaling limit the local time distribution P\oc{T\t) tends to the 
following asymptotic form 



Pi.c{T\t) ~ exp 



(40) 



where $(r) is a large deviation function. 

To compute the large deviation function we first substitute 
this presumed asymptotic form of P\oc{T\t) given by Eq. ( I40t 
in the Laplace transform G(a) ~ /g e~"*Pioc(r|t) and 
then make a change of variable r = T/t m the integra- 
tion. The resulting integral can be evaluated in the large T 
limit by the method of steepest descent, which gives G{a) ~ 
exp[-TW(a)] where W{a) = minr[{a + <I>(r)}/r]. Com- 
paring this result with Eq. ( I28> gives 



$(r) 



A(a), 



(41) 



where A(a) is given by Eq. \26l . Thus \{a) is just the Leg- 
endre transform of $(r). Inversion of this transform gives the 
exact large deviation function 



$(r) — maxQ,[— a + r\{a)], 



(42) 



with \{a) given by Eq. J26> . This is a general result valid for 
any confining potential U {x). 

We will now explicitly compute the large deviation $(r) for 
the particular potential given by Eq. (|5} with /i < and cr = 0. 
Substituting the corresponding force F{x) = — sign(a;) in 
Eq. (I23> and solving the resulting differential equations with 
the boundary conditions y+{x 00) and y-{x 
—00) ^ we get 



2/±(a;) = 2/±(0)exp T + Va^^ + 2a ) a; 



(43) 



Substituting these results in Eq. i26\ we get A(a) = — + 
2a. From Eq. (I42t one then gets the large deviation 



function 



(44) 
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It turns out that for this particular form of the force F{x) = 
— sign(a;), the Laplace transform in Eq. j28> can be in- 
verted to get the local time distribution Pioc{T\t) exactly for 
all T and t. The calculations are presented in appendix 1X1 
We find that in the large t limit, the distribution reduces to the 
asymptotic form 



: exp 



t /T 

2 It 



(45) 



near the mean (T) = which verifies the result obtained 
above by the large deviation function calculation. 

In fact, the limiting Gaussian form of the distribution of the 
local time near its mean value is quite generic for any stable 
potentials (where the system eventually becomes ergodic) and 
is just the manifestation of the central limit theorem. From the 
definition, T - (T) = J*{S[x{t')] - {S[x{f)])} df , it follows 
that when T ^ (T), the random vmiahles S[x{t')]- {S[x{t')]) 
at different times t' become only weakly correlated. Then in 
the limit when t is much larger than the correlation time be- 
tween these variables, one expects the central limit theorem 
to hold which predicts a Gaussian form for T near its mean 
value (T). 



where K^{x) is the modified Bessel function of order v ll44ll 
and K^^{x) — K^{x). Averaging Eq. (I28> over disorder we 
finally get the exact formula 

However, it is not an easy task to invert the Laplace transform 
to get the exact distribution P\ociT\t) for all T and t. In the 
following subsections we will extract the asymptotic behav- 
iors of Pioc{T\t), for the three cases, when the deterministic 
part of the potential is: (i) flat corresponding to /i = 0, (ii) un- 
stable corresponding to /x > and (iii) stable corresponding 
to /X < 0. 



A. Flat potential {/i — 0) - Sinai model 

We first consider a particle diffusing in the continuous Sinai 
potential, i.e., /i = in Eq. Q- Our aim is to find out how 
this random potential modifies the behavior of the local time. 
In this case substituting q{T) and q{0) from Eq. ( l48l with 
/i = in Eq. i49i we get the Laplace transform of the disorder 
averaged local time distribution as 



IV. LOCAL TIME WITH DISORDER (a > 0) 

So far we have considered the case where the random part 
of the potential was not present. In this section we will study 
the effect of the randomness by adding a random part to the 
potential. In particular, we will consider the diffusive motion 
of the particle when the force F{x) is given by Eq. (|6} with 
cr > 0. 

Equation ( I28t still remains valid for each realization of the 
force F{x), i.e., for each realization of {^(a;)}. Our aim is 
to compute the average of the pdf of the local time Pioc{T\t) 
over the noise history {(,{x)}. From Eq. ( I28> , one needs to 
know the distribution of X{a) = [^-(0) — z-^-{0)]/2, which is 
now a random variable since F{x) is random. The variables 
— z+(0) and z^{Q) are independent of each other and there- 
fore their joint probability distribution factorizes to the indi- 
vidual distributions. The calculations of these distributions 
are presented in appendixiBl Using the distributions of z±{0) 
from Eqs. ( IB9> and JBl 1> respectively with a± — a, one gets 



exp[-A(a)T] = 



QjT) 

m 



(46) 



with 



q{T) = / tu^/'^^^cxp 

^0 



If, , 2a 

2(7 w 



^2a{l + aT) 



dte-'''Pi,,{T\t) 



' ^2a{l + aT) 



(50) 



We will now consider the interesting limit where both t and T 
are large, but the ratio y = T/tis kept fixed. This corresponds 
to taking the limit a ^ with aT = s keeping fixed. In this 
limit 



2a 



1 



■ logo; 



(51) 



and 



(52) 



Therefore substituting t = s/ay and T = s/ain Eq. (I50> . in 
the limit a ^ we get 



dye' 



s ( s 






-) 




ay J 



J2 a„-^0 



2s 



(53) 



log ads " \\/a 

The above equation suggests that, in the limit t oo and 
T oo, while their ratio T/t is kept fixed, P\oc{T\t) should 
have the scahng form 



Pxoc{T\t) 



(54) 
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Now substituting the above form in Eq. J53t and making the 
change of variable y = 1/y, we obtain after straightforward 
simpHfication 



dye' 



-sy 



him 



y2 



(55) 



Note that the right hand side of the above equation is simply 
the Laplace transform of the function Therefore 
by using the identity 



oo g— a"/4ii) 



' duj = 2Kq (ay/s) 



(56) 



and the convolution property of Laplace transform, we can in- 
vert the Laplace transform in Eq. i55\ with respect to s. Invert- 
ing the Laplace transform and after simplification we finally 
get 



/i(l/y) = 2y 



1/2 



dx 



x{l~x) 



exp 



1 



2ayx{l — x) 



Therefore the scaling function /i (y) is simply given by 



fiiy) 



1/2 



dx 



y Jo x{l - x) 



exp 



y 



2ax{l — x) 



By substituting x(l — x) — 1/z gives 



/i(y) 



dz 



y Ji y/z{z - 4) 



exp 



(57) 



(58) 



(59) 



where the integral can be evaluated exactly which finally 
gives the scaling function in Eq. ( I54t as 



h{y) = -e-y^'^Koiy/a). 

y 



(60) 



However the scaling given by Eq. (I54> breaks down for very 
small y (very small T) when y <^ a. The scaling function is 
displayed in Fig.|3] In the large y limit, using the asymptotic 
the behavior K,, (x) ~ yj-nj^x e^^ from Eq. ( l60l we find that 

h{y)^V2^y-y'e-^y/-. 



B. Unstable potential {fi > 0) 

In this case the behavior in the limit t ^ oo can be obtained 
by either setting a = in the integral in Eq. (I47> or taking the 
a ^ limit in K,y{.) in Eq. ( I48t . which gives 



g(T)^r(M/a)(2a)'^>(l + ar)-^/^ 



(61) 



where r(a;) is the Gamma function |45]. Substituting q{T) 
and q{Q) in Eq. ( I49> and inverting the Laplace transform with 
respect to a gives 



Pioc(T|t) = 2^i{l + aT) 



-2/j/cr-l 



(62) 




FIG. 3: The scaling function fi{y) in Eq. <54t . The solid line is 
plotted by using Eq. <60> . and the dash line is plotted by using the 
limiting form /i (y) 



/2nay-^/'^e''^y^'' as y - 



i.e., in the limit t oo, the distribution Pioc{T\t) tends to 
a steady state distribution Pioc{T) for all T > 0. The dis- 
order averaged local time distribution has a broad power law 
distribution even though for each sample the local time has 
a narrow exponential distribution (see Ea.(l33> in Sec. iniB> . 
This indicates wide sample to sample fluctuations and lack of 
self-averaging. 



C. Stable potential {fi < 0) 

In this case substituting q{T) and q{0) from Eq. (I48> in 
Eq. (I49> and denoting u — \ii\/a we get 



1 



dT 



(63) 



We consider the scaling limit where both t and T are large, 
but their ratio y = T/t is, kept fixed. This corresponds to 
taking the limit a ^ with keeping aT = s fixed, which 
gives the following limiting forms 



a 



\/2a\ T(v) ( cr\/2 




(64) 
(65) 

(66) 
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Substituting the above limits on the right hand side of Eq. (16 3> 
and making change of variables t = s /ay and T = s/a on 
the left hand side, it is straightforward to get 



dye~ 





f s 








-) 


ay' 




ay J 



d 



/2s 



(67) 



{2aYV^{u) ds 

in the limit a 0. This suggests the limiting form for 

Ploc(T|t) 



1 



Pxoc{T\t) ^ -/2(rA), 



(68) 



in the scaling limit t ^ oo and T ^ oo with a fixed ratio y = 
T/t. To compute the scaling function we substitute the above 
scaling form in Eq. ( I67t . and make the change of variable y = 
1/y. Then Eq. (|^} simplifies to the Laplace transform 



dye' 



f2{l/y) 



y2 



/2s 



, (69) 



which can be inverted with respect to s, by using the identity 



(?) '-^e—d^^2s^/^K4aV-s) 



(70) 



and the convolution property of the Laplace transform. After 
simplification, the inverse Laplace transform gives 



/2(y) 



1/2 



dx 



y 



X exp 



2(Tx(l — x) 



(71) 



By making a change of variable x{l — x) — 1/z in the above 
integral, it can be presented in the form 



/2(y) = 



2y2 



(z-4) 



-1/2 



X 



exp (-^^) dz, (72) 



which now can be expressed in more elegant forms 14511 as 

2iy-l 



f2{y) 



' 2^ ■ 









(J) e-2^/^{/(l/2,l + ^.,2j//a), 



(73) 

where [/ (a, 6, a;) is the Confluent Hypergeometric Function 
of the Second Kind (also known as Kummer's function of the 
second kind) which has the following limiting behaviors: 



[/(l/2,l + zy,a;) « ^^x"'' forsmafla;, 



U(\l2,\^v,x)^ ^ 

V X 



for large x. 



(74a) 
(74b) 



0.5r 



0.3- 
b 0.2 - 



\v > 1 



W < 1 



1 2 3 4 5 

y/(^ 

FIG. 4: The scaling function /2(j/) in Eq. <68> . plotted by using 
Eq. <73> . !y = ImI/''"- 



The scaling function /2(y) is displayed in Fig.|4] Using the 
limiting behaviors from Eq. J74> . one finds that the scaling 
function decays as /2(y) ~ y(4''-3)/2g-2y/(T for large?/. For 
small y, the scaling function behaves as f2{y) ^ y'^^^, which 
increases with y for ly > 1, however, diverges when y 
for 1/ < 1, a behavior qualitatively similar to the Sinai case 
(see Fig. O. For i' < 1, the disorder wins over the strength 
of the stable potential. In that situation when the particle gets 
trapped in the wells of the random potential, the weak external 
deterministic force often can not lift it out of the well and 
send towards the origin. Therefore, the scaling function f2{y) 
carries very large weight near y — (which corresponds to 
very small local time T for a given observation time t). 

Note that, for the particular value ly — 1/2, the scaling func- 
tion has a simple form /2(y) = \/2/7ray exp{—2y/a). 



V. INVERSE LOCAL TIME WITHOUT DISORDER (a = 0) 

The inverse local time means how long one has to observe 
the particle until the total time spent in the infinitesimal neigh- 
borhood of the origin is T. The double Laplace transform of 
the pdf of the inverse local time is obtained by simply putting 
X = in Eq. il9i . Corresponding u(0) in Eq. il9i . which is 
nothing but the double Laplace transform of the pdf of local 
time, has already been evaluated in Sec. Hu] and is given by 
Eq. (|23. Substituting w(0) and replacing /(t|r, 0) with the 
pdf of the inverse local time Iioc{t\T), after straightforward 
simplification, for x = Eq. ( I19t reads 



poo poo 

/ dte-"' / dT e-P^ Iioc{t\T) = 
Jq Jq 



1 



p + X{a) ' 



(75) 



where A(a) is given by Eq. ( I26> , which depends on the force 
F{x) through Eq. ( I23t . Inverting the Laplace transform with 



10 




t/T 



FIG. 5: The pdfs of the inverse local time for stable (/i = —1/2), flat 
(/^ = 0) and unstable {^l = 1 /2) potentials, plotted using Eq. |7i 
and r = 2. 



respect to p gives the general formula 

POO 

/ e-"*/ioe(t|r) = exp [-\{a)T] , (76) 

valid for arbitrar y fo rce F{x), a result known in the mathe- 
matics literature I17ll4al . 

We first consider the pure case where the force given by 
Eq. (|6} with a ~ Q. Substituting solutions of Eq. ( I23> 
for F{x) — /isign(j:) in Eq. M6\ we obtain \{a) = /i + 
+ 2a. Now using this \{a) in Eq. M6\ and making a 
change of the parameter a = /3 — /^^/2 we get 



it decays exponentially /ioc(i|T) exp(— /x^</2). On the 
other hand when the potential is unstable, /i > 0, as we see 
from Eq. M9\ . the distribution /ioc(^|T') is not normalized to 
unity. In this case the particle escapes to ±oo with probability 
(1 — e~^^^) and Eq. (I78> gives the distribution only for those 
events where the particle does not escape to ±cxd. Therefore 
for /X > 0, it is appropriate to represent the full normalized 
distribution as 



: exp 



2t 

[l-e-^''^)S{t-oo). (80) 



Note that the second term does not show up in the Laplace 
transform of Iioc{t\T) with respect to t. 



VI. INVERSE LOCAL TIME WITH DISORDER (a > 0) 

In this section, we switch on the disorder by considering 
CT > in the force given by Eq. (|6j. In the presence of disor- 
der, taking the disorder average ofEq. ^ gives 



dtf 



*/,oe(i|T) = exp [-A(a)r], 



(81) 



withA(a) = [z_(0)-z+(0)]/2, where -z+(0) and z_(0) are 
independent random variables, whose distributions are given 
by Eqs. ( IB9> and ( IBl U respectively with a± = a. The object 
exp [— A(a)T] on the right hand side of Eq. J8U . has already 
been evaluated in Sec. IIVI which is given by Eq. i46i . In 
the following subsections we will determine the behavior of 
Iiocit\T) in the scaling limit t oo, T oo, while keeping 
their ratio x ^ t/T fixed, for the three qualitatively different 
cases: (i) /i = 0, (ii) /i > and (iii) ii < 0. 



dtf 



-pt 



^y'^^Iioc{t\T) 



(77) 



where the right hand side is simply the Laplace transform of 
*-/'^Iy^{^t\T) with respect to t. The Laplace transform can 
be inverted to obtain the exact pdf of the inverse local time 



hoc{t\T) 



T 



exp 



%/2^ 

with the normalization condition 



2t 



(78) 



A. Flat potential (/i = 0) - Sinai model 



Following the analysis of Sec. II V Al in the limit a 
with keeping oT — s fixed. 



exp [-\{a)T] 



log a 



2 -"^0 



'2s 



(82) 



Therefore, substituting t = xs/a and T = s/a, in the limit 
a — > 0, Eq. ( 18 1> reads 



/,oc(t|r)dt-e-(^+i^i)^ 



1 for ^ < 0, 
g-2pT fQj. ^ > 0, 

(79) 

which is simply obtained by putting a = in Eq. ( I76> . As we 
infer from Eq. (I78> . although in the limit t ^ the inverse 
local time distribution Iioc{t\T) ~ exp(— r^/2t) is indepen- 
dent of /i, for large t it depends on the nature of the poten- 
tial, as shown in Fig |5] While in the absence of any force, 
i.e., ^ — the inverse local time distribution has a power-law 
tail IiociAT) ^ for the stable potential, i.e., ^ < 0, 



dx f 





S 


/ sx 




— I\oc 




-) 




a 


V a 


a / 



This suggest the scaling form 



/ioc(i|r) 



1 



log a 



(83) 



(84) 



in the limit t — > cx3, T ^ cx) but keeping x — t/T fixed. Sub- 
stituting this scaling form in Eq. ( I83> . after straightforward 
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b 




FIG. 6: The scaling function gi{y) in Eq. <84t . The solid line is 
plotted by using Eq. <86> . and the dash line is plotted by using the 
limiting forms: gi{x) ~ \/2nax^^^^ exp(— 2/ax) for small x and 
gi(x) ~ 2 log(crx)/x for large a;. 



simplification one obtains 



2s 



(85) 



Now direct comparison of the above equation with Eq. (I55> 
gives gi{x) — fi{l/x) jx^, where ]\{x) is given by Eq. ( l6Qt . 
Substituting /i(l/x) one obtains the scaling function 51 (a;) 
as 



9\ 



(x) = - e-'^-^Ko (l/ax) , 



(86) 



which is displayed in Fig.|6l The scaling function increases as 
gi{x) ~ \/2tt<7x^^/^ exp{—2/ax) for small x and decays as 
gi (x) ^ 2 \og{ax) /x at large x. 



B. Unstable potential (fj, > 0) 



In this case the right hand side of Eq. ( I8U is given by 

exp[-A(a)r] = il+aT)-" V2 ^ /^/ A ^' ^^^^ 

with 1/ = n/a. Putting a = in the above equation gives 
the normalization condition, Iiocit\T) dt = {1 + aT)^^'^, 
which implies that for the unstable potential, where the force 
is repulsive from the origin, the particle escapes to ±00 with 
probability 1 — (1 + aT)^'^'^, and the disorder averaged 
pdf I\oc{t\T) obtained by inverting the Laplace transform in 
Eq. ( I8II 1 represents only those events where the particle does 
not escape to ±00. 



Now in the limit of a — > with aT = s keeping fixed, one 
gets 



2a 



r>) y./2 ( ^ 



Therefore Eq. j87l i becomes 



'2s 



(88) 
(89) 



exp [-\{a)T] 



{2a'^YT^{v) 



'2s 



(90) 



In the corresponding limit T 00, t 00, but keeping 
their ratio x = t/T fixed, using the scaling form 



1 



hoc{t\T) = ,^^g2{t/T) 
in Eq. J81> one finally arrives at the Laplace transform 



(91) 



e '^^g2{x)dx 



{2a^YV^{v) 



'2s 



(92) 

The Laplace transform can be inverted with respect to s to ob- 
tain the scaling function 52(2;) and in fact the inversion has 
akeady been done in Sec. IIV CI Comparing the above equa- 
tion with Eq. ( I69> readily gives .92(2^) — (^^"^^ f2{^ / x) / x"^ 
where f2{x) is given by Eq. ( I73> . Substituting /2(l/x) gives 



52(2;) 



-2/<jx 



Uil/2.l + ,y,2/ax), (93) 

[axj'^'^^^ 



where U{a,b,x) is the Confluent Hypergeometric Function 
of the Second Kind, whose its small and large x behaviors 
are given in Eq. i74\ . The scaling function g2{x) is dis- 
played in Fig. The scaling function increases as g2{x) ^ 
exp(— 2/(Tx) for small x and eventually decreases for large x 
as g2{x) ~ l/x'^" . In particular, for 1/ — 1/2 it has a very 
simple form 52 (a;) = \/2/7ra^x'^^/^ exp{~2/ax). 



C. Stable potential {fi < 0) 



Following the analysis of Sec. IIV CI in the limit a ^ 0, 
keeping aT — s fixed one gets 



exp [-X{a)T] = 



(2(7)^2(1/) 



s'^^'K. 



'2s 



, (94) 



with ly — ImI/o"- 

On the other hand, in the corresponding limit T ^ oo,t 
00, but keeping t/T — x fixed, using the scaling form 



1 



Iioc{t\T) = -g3{t/T), 



(95) 
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FIG. 7: The scaling functions 32(2) in Eq. <91> plotted by using 
Eq.^.iy = \^i\/a. 



one gets 



/>oo />oo 

/ dte-°'*Iioo{t\T) = / dx( 
^0 ^0 



■^53(2;), (96) 



with s = aT. Therefore, in this scaling limit Eq. ( 18 U be- 
comes 



rf-^-^^53(x) = ^2a)'^r^M 



/2s 



H 2 



Now comparing the above equation with Eq. J92> one gets 



(97) 



53(2;) 



2^ 



(98) 



where the scaling function 32 (2;) is given by Eq. ( I93> and dis- 
played in Fig.0 While Iioc{t\T) has the same scaling function 
(up to a multiplicative factor of a^'^) for both stable and unsta- 
ble potential, the physical behaviors, however, are quite differ- 
ent. For the stable potential, Iiocit\T) is normalized to unity. 
Note that the scaling function 53(2;) becomes narrower as one 
increases as expected since the particle becomes more lo- 
calized near the origin. For the unstable potential, on the other 
hand, the weight of Iioc{t\T) decreases as {aT)~'^'^, as one in- 
creases ly, as expected since when the repulsive force from the 
origin becomes stronger, the particle escapes to ±00 with a 
higher probability. 



VII. OCCUPATION TIME WITHOUT DISORDER (a = 0) 

In this case V{x) = 0{x), corresponds to T in Eq. (|2jl be- 
ing the occupation time in the region a; > and P{T\t, 0) in 
Eq. being Poc<:[T\t) — the pdf of the occupation time for 
a given observation time window of size t and the initial po- 
sition of the particle x(0) = 0. Again as before, we need to 



solve the differential equation ( I14> for x > and a- < sep- 
arately and then match the solutions at x = 0. The matching 
condition for the slope of the solutions is obtained by integrat- 
ing Eq. (I14> across a; = 0. Thus the matching conditions are 

u+(0) = u_(0) = u(0), and u\{Q) = m'_(0), (99) 
where u±{x) satisfy the following differential equations: 

]^ul{x) + F{x)u'^{x) - {a+p)u+{x) = -1, (100) 

for a; > 0, and 
1 



u"_ {x) + F{x)u'_ {x) — au^ {x) 



(101) 



for a; < 0. The boundary conditions of u± {x) when x — > ±cx) 
are obtained from the fact that if the starting position goes 
to ±00, the particle will never cross the origin in finite time: 

P{T\t, a; -> 00) (5(t - T) and P{T\t, x -00) = (5(T), 
and hence from Eq. M2>\ 



M+(oo) = , and 7i_(— 00) ~ —. 

a + p a 



(102) 



Writing u+{x) — \/{a + p) + B^y^{x) and u^{x) = 
1/a + B^y^{x), we obtain the homogeneous differential 
equations for y± (x) as 



ly'lix) + F{x)y'+{x) ~ {a+p)y+{x) = 0, 



for a; > 0, and 
1 



y'L[x)+F{x)y'^{x)-ay^[x)^Q, 



(103) 



(104) 



for a: < 0, with the boundary conditions (.x ^ 00) = and 
j/_ (a; — > -co) ~ 0. The constants B± are determined by the 
matching conditions given in Eq. (|99}, which can be rewritten 

as 

^— + B+y+(0) = - + B_y_(0) =m(0), (105a) 
a + p a 



B+y'+{0) = B-yUO). 



(105b) 



Eliminating the constants from Eq. (I105> . we obtain the 
double Laplace transform of the pdf of the occupation time 



Jo Jo 

h{a,p) i2{a,p) 



(106) 



a+ p 



where 



z_(0) 



z_(0)-z+(0) 
-z+(0) 



Lz_(0)-z+(0) 



(107) 
(108) 
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and z±{x) = y±{x)/y±{x). Note that 

ii{a,p) + £2{a,p) 1. 



(109) 



Putting p ~ 0, in Eq. ( I106t gives u(0) = 1/a, and hence in- 
verting the Laplace transform with respect to a readily verifies 
the normalization 

/ Pocc(T|<)dT= 1. (110) 

For any symmetric deterministic potential the distribution 
of the occupation time is symmetric about its mean (T) = t/2, 
i.e., Pocc{T\t) = Pocc{t - T\t). Then, it follows from this 
symmetry that 



£i(a+j3, -p) ^£2{a,p). 



(Ill) 



In other words, the double-integral in Eq. ( I106t remains in- 
variant under the following simultaneous replacements: (a + 
p) — *■ a and a {a + p). Thus under these replace- 
ments one must have z+{Q) — > — z_(0) and vice versa, which 
also implies that z+(0) = — z_(0) for p — 0. Equivalently, 
£i(a, 0) = £2 (a, 0) = 1/2, which also directly follows from 
Eqs. ( fTU9l l and (fTTTl . 

Therefore if one splits the distribution function into two 
parts: Pocc(T|i) = RL{T\t) + RR{T\t) such that 



Jo Jo 

/ dte-"* / dT e-P'^ RR{T\t) = 
Jo Jo 



a 

a+p 



(112) 
(113) 



then it follows from the above discussion that RL{t — T\t) = 
Rii{T\t). This symmetry of the distribution will come handy 
later Moreover, putting p — Q and inverting the Laplace trans- 
forms with respect to a gives the normalization for each part 
separately 



/ RL{T\t)dT= f RR{T\t)dT=]-. (1 
Jo Jo ^ 



14) 



As an example, we first consider the pure case: cr = in the 
force given by Eq. (|6j. For F{x) = /^sign(x), the solutions 
of Eqs. ( fT03t and(fT04l are obtained as 



, (115) 



= 2/+(0)exp - + V/^^ + 2(a + p)^ x 
for x > 0, and, 

[( 



V~{x) = y-(0)cxp 



2a] X 



(116) 



for a; < 0. These give the expressions for z±{0) 
yi(0)/j/±(0) as 



z+(0) = - + 2(a+p) 



z_(0) = pi+vJj^+2c 



(117) 



In the following subsections we will consider the three differ- 
ent cases: (i) /i = 0, (ii) /i > and (iii) /i < 0. 



A. Flat potential {fi — 0) 



For /i = 0, using z^{0) — — \/2{a + p) and z_ (0) — v2a 
from Eq. ( I106> we get 



/ die-"* / dre-p^Foec(r|0 ^ , , (1 

Jo Jo y^a{a+p) 



18) 



Inverting the double Laplace transform with respect to p and 
then with respect to a finally reproduces the well known 
Levy's "arcsine" law (III] for the pdf of the occupation time 
of an ordinary Brownian motion. 



Pocc(r|i) 



<T <t. 



(119) 



The distribution Pocc{T\t) diverges on both ends T ~ Q and 
T ~ t, which indicates that the Brownian particle "tends" to 
stay on one side of the origin. 



B. Unstable potential (/i > 0) 

Since for fi > 0, the force is repulsive from the origin 
X — 0, one would expect the occupation time distribution 
to be convex (concave upward), with minimum at T = t/2. 
Now in the limit of large window size t ^ 00, the part of 
the distribution Pocc(T|<) to the left of the midpoint T = t/2 
approaches i?i(r|t), as the midpoint itself goes to 00. 

By making a change of variable z — at, it follows from 
Eq. (tmi 



dz e 



dTe~P^RL{T\z/a) = li{a,p). (120) 



Now the large t limit of RL{T\t) can be obtained by taking 
a — > in the above equation, where one realizes that i?i(r|t) 
approaches a steady {t independent) distribution, i?L(T|< 
Qo) — > Rl (r), whose Laplace transform is given by 



dTe'P^RL{T) = h{Q,p), 



(121) 



where ii{Q,p) is obtained from Eq. ( I107> . by using 2:±(0) 
from Eq. ( II 17L which gives 



^1(0, p) = 



2/i 



3Ai+ Va^^ + 2p' 



(122) 



The above Laplace transform can be inverted with respect to 
p, which gives 



Rl{T) = ^l^/2e- 



3/i 



1 



V2 V 2 ) V\/2 
with the normahzation /J^ Rl{T) dT = £i{0, 0) = 1/2. 



T erfc 



3^ 



(123) 
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The limiting behavior of the distribution is given by 



RUT) 



(124) 



for small T and decays exponentially for large T, 



Rl{T) 



V2 



-l?Tj2 



(125) 



C. Stable potential {p. < 0) 



As we discussed earlier in Sec. 1111 CI in the context of 
the local time, for generic stable potential U{x) the system 
eventually becomes ergodic at large t and hence the average 
{6[x{t)]) approaches its stationary value {9[x{t)]) 
where Z ~ g-2C/(x) is the equilibrium partition func- 
tion and Z+ ~ g-'^u(x) jj^g restricted partition 
function. Therefore, for large t the average occupation time 
(T) = /o (6'[x(t')]) dt' scales Hnearly with t 



(126) 



Note that when the potential U{x) is symmetric about zero, 
the average occupation time (T) — t/2 for all t. 

From the definition, T - (T) ^ Jo{d[x{t')] - 
{e[x{f)])}dt', it follows that when T (T>, the random 
variables 9[x{t')] — {6[x{t')]) at different times t' become only 
weekly correlated. Then in the limit when t is much larger 
than the correlation time between these variables, one expects 
the central limit theorem to hold, which predicts a Gaussian 
form for the distribution of the occupation time T near the 
mean value (T), 



Pocc(r|i) -exp 



where the variance cr^ = (T^) — (r)^ can be obtained from 
the Laplace transform of the moments 



{T')e-°'Ut= (-1) 



,a"u(o) 



(128) 



p=0 



with u{0) given by Eq. jl06> . 

For the particular attractive force F{x) ~ — sign(a;), us- 
ing z± (0) from Eq. il 17> in Eq. (I106> and taking derivatives 
with respect to p in we get 



du{0) 



dp 



dp^ 



p=0 



p=0 



1 
1 

2^ 



O 



(129) 



(130) 



Therefore inverting the Laplace transform in Eq. ( I128> with 
respect to a immediately gives (T) = t/2 for all t, and 



VIII. OCCUPATION TIME WITH DISORDER (a > 0) 

Now we consider the occupation time when the disorder is 
switched on: > in Eq. (|6}. Our aim is to calculate the 
disorder averaged Pocc{T\t). As one realizes from Eqs. ( I106t . 
( fTOTl l and ( fTOSt . to calculate 'P^^f[t) one needs the distri- 



bution of — z+(0) and z_(0), which are given by Eqs. (IB9> 
and ( IBl 1> with a+ = a + p and a_ = a respectively. In 
the following subsections, we will consider the three cases; 
(i) /i = 0, (ii) ^ > and (iii) fi < 0. 



A. Flat potential {fj, = 0) - Sinai model 

We first consider the diffusive motion of a particle in a con- 
tinuous Sinai potential, where the potential itself is a Brown- 
ian motion in space. In the limit of large window size t the 
left half of the disorder averaged pdf of the occupation time 
i?L(T|i) for < T < i/2 is obtained by taking the disor- 
der average in Eq. il 12> . The right half of the distribution 
for t/2 < T < t is just the symmetric reflection of the left 
part. The detailed calculations for Ri^{T\t) are presented in 
appendixicl 



We find that i?i(T|i) has a large t behavior 

1 



logt 



Rin 



(131) 



where the function R{T) is independent of t. The limiting 
behaviors of R{T) are given by 



R{T) 



V2(j 

7^' 



as T ^ and 



(127) for large T. 



R(T) , 



B. Unstable potential {/i > 0) 



(132) 



(133) 



For /i > 0, we find that disorder does not change the asymp- 
totic behavior of the distribution for the pure case qualita- 
tively. The calculations are presented in appendix|Dl We find 
that in the limit t ^ oo the left half of the disorder averaged 
occupation time distribution tends to a t independent form 



RL{T\t) ^ Rl{T). 



(134) 



In fact the small T limit of Rl{T) remains same as in the pure 
case 



Rl{T) 



mV2 



(135) 



For large T, the distribution Rl{T) stiU decays exponentially 



Rl{T) 



-bT 



(136) 
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where the decay coefficient b is, however, different from the 
pure case (see Eq. lD17> . 



C. Stable potential (fi < 0) 

This particular situation, where one finds the interplay be- 
tween two competing processes, is a very interesting one. On 
one hand, as we discussed in Sec. I VII Cl the stable potential 
in the absence of the disordered potential makes the system 
ergodic in the large t limit, and as a result the pdf of the oc- 
cupation time is peaked at = t/2 and decays fast away from 
it. On the other hand, as we discussed in Sec. IVlll Al without 
any underlying deterministic potential the disorder averaged 
pdf of the occupation time is convex (concave upward) with a 
minimum at T = t/2 and diverges at the both ends T ^ and 
T ^ t. Therefore, if both the stable potential and disordered 
potential are included, as their relative strength i/ = |/i|/cr is 
varied, one expects a phase transition at some critical value i^c 
where the system looses ergodicity. 

In the scaling limit where both t ^ oo and T ^ oo, but 
their ratio y — T/t is kept fixed, we find that the disorder 
averaged pdf of the occupation time has a scaling form 



1 



Pocc(T|t) = -fo{T/t) 



(137) 



The calculation of the scaling function fo{y) is presented in 
appendix|EJ where we find the Beta law 



foiy) 



B{iy, I/) 



0<y<l, (138) 



where v = \^\/cr and B{v,v) is the Beta function 13 • Now 
if one tunes the parameter v by varying either /i or the disorder 
strength cr, the distribution Pocc{T\t) exhibits a phase transi- 
tion in the ergodicity of the particle position at = 1 (Fig ID- 
For V < Vc the distribution fo{y) in Eq. (I137> is convex with 
a minimum at y ~ 1/2 and diverges at the two ends y — 0,1. 
This means that particle tends to stay on one side of the origin 
such that T is close to either or t. In other words the paths 
with small number of zero crossings carry more weight than 
the ones that cross many times. For v > Vc the scenario is 
exactly opposite, where fo (y) is maximum at the mean value 
y = 1/2 indicating that particle tends to spend equal times 
on both sides of the origin x = 0, such that paths with large 
number of zero crossings, for which T is closer to t/2 carry 
larger weight. Similar phase transition in the ergodicity prop- 
erties of a stochastic process as one changes a parameter, was 
first noted in the context of diffusion equation 1 19|, and later 
found for a class of Gaussian Markov processes |20] and in 
simple models of coarsening f47','48'|. 

A very interesting observation about Eq. ( I138> is that for 
1/ — 1/2, the result is same as Levy's result for the one- 
dimensional Brownian motion given by Eq. ( II 19> . It seems 
as if the attractive force cancels the effect of disorder exactly 
at 1/ = 1/2. However, this is no more true in the context of 
the local time. 




FIG. 8: The scaling functions fo{x) in Eq. <137t plotted by using 
Eq. (n^ . 



IX. INVERSE OCCUPATION TIME WITHOUT 
DISORDER (cr = 0) 

In this case I{t\T, 0) in Eq. ([19) is replaced with Iocc{t\T), 
which is the distribution of the time t needed to observe the 
particle with a starting position a; = 0, until the total amount 
of time spent on the positive side x > is T. Correspond- 
ing u(0) in Eq. ( I19> with x = 0, which is the double Laplace 
transform of the pdf of the occupation time, has akeady been 
evaluated in Sec. IVlll which is given by Eq. ( I106> . Substitut- 
ing u(0) in Eq. il9i gives 



.m = '-^, (139) 

a+p 

where £2(0:, p) is given by Eq. ( I108> . Comparing the above 
equation with Eq. dl 13> . one can infer that Iocc{t\T) and 
Rii{T\t) have the same functional form, i.e., Ioccit\T) = 
Rii{T\t) and especially for symmetric deterministic potential 
Ioccit\T)=RRiT\t) = RL{t~T\t). 

It is useful to present the above equation in the following 
form. 



dz ( 



dr ( 



/3 



(140) 

which has been obtained by substituting p — (3 — a in 
Eq. ( I139> and subsequently making the change of variables 
f3T — z and t ^ t — T. On the right hand side, we have 
substituted £2 (a, f3 — a) = £i{P,a — j3), using Eq. (fTTTl and 
£i (a, p) is given by Eq. ( I107t . Now by taking the limit /3 ^ 
in Eq. ( I140> . one obtains the large T behavior of /occ(^|r). 

For the pure case, cr = in Eq. (|6}, we have already 
obtained 2;-|-(0) in Sec. IVlll which are given by Eq. dl 17l i 
and hence we can evaluate £i{a^p) and l2{a^p) by using 
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FIG. 9: The pdf of the inverse occupation time for simple Brownian 
motion. 



Eqs. ( llQ7t and ( I108> respectively. In the following subsec- 
tions we will analyze the behavior of Iocc{t\T) for the cases 
(i) ^ = 0, (ii) /i > and (iii) /i < 0. 



A. Flat potential (/i = 0) 

For fi = 0, which is the case of a simple Brownian motion, 
z+(0) = —y^2{a + p) and z^{0) = \f2a. Therefore using 
Eq. ( fTO . from Eq. ( fll^ we get 



dTe"P^ / dte-"*Iocc(.t\T) 



1 



(141) 



Now inverting the Laplace transform with respect to p gives 

/oo 
die-"*/occ(i|r) = erfc(V^), (142) 

and further inverting the Laplace transform with respect to a 

gives 



Iocc{t\T) = 



Vt 



9{t~T). 



(143) 



with the normalization condition Iocc{t\T) dt — 1, which 
is readily checked by putting a = in Eq. (I142> . The in- 
verse occupation time has non-zero support only for t > T,as 
shown in Fig.|9] 

Note that, since RR{T\t) = Ioccit\T) = VT/TTtVt^ 
and B.L{T\t) = Rn{t - T\t) = Vt^^/TrtVT, adding 
the t wo parts, Pocc(T|t) = RL{T\t) + RR{T\t) = 
l/7r^r(t - T), one recovers Eq. ( ITT91 . 



B. Unstable potential {/i > 0) 

For /i > 0, although one can invert the Laplace transform 
in Eq. ( I139> with respect to p exactly, the other Laplace trans- 
form with respect to a can be inverted only in the large T 



limit. Therefore to keep the presentation simpler, we will con- 
sider the large T behavior of /occ(i by analyzing Eq. jl40t 
in the limit (3^0. 

By using z± (0) from Eq. il 17> in Eq. ( llQ7t . one gets 



^i(0,a) 



2^ 



+ 2a 



(144) 



Therefore, Eq. J 1401 1 suggest that Ioccit\T) should only de- 
pend on the difference {t — T) at large T 

Ioccit\T)=hit-T). (145) 

Substituting this form in Eq. (ll4Qt in the limit /? — > gives 



■/i(t) =fi(0,a), 



(146) 



where putting a = gives the normalization Ii (r) dr = 
ii{0, 0) = 1/2, indicating that the particle can escape to — oo 
with probability 1/2 for the unstable potential (the force is re- 
pulsive from the origin). Now inverting the Laplace transform 
with respect to a gives 



/i(r) =/iV2e-^'^/2 



1 



3/^ /V A f ^^^^ r 



(147) 



(148) 



for small r = [t — T) and decays exponentially for large 

T={t^T), 



The limiting behavior of this distribution is given by 
Ii{t) 



hir) 



^/2 e-^"-/^ 
9m\AF r3/2 • 



(149) 



C. stable potential (/i < 0) 

It is reasonable to consider the difference variable t — T 
instead of t, as t > T. Therefore, we write 

Ioccit\T)=l2it-T,T). (150) 

Substituting this form and p = (3 ~ ain Eq. ( 1 1391 1 one gets 



dTe 



dre-""/2(T,T) 



ei{P,a-f3) 



(151) 

where we have substituted £2(0;, f3 — a)= ii{(3, a — (3) on the 
right hand side, using Eq. (^Hj. Using ^±(0) from Eq. il 17> 
for /i < 0, in Eq. ( I107t gives 



vV + 2/3+ Vm^ + 2^-2|^| 



(152) 
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Therefore, taking the small /3 limit in Eq. ( I151> gives 



, (153) 



and inverting the Laplace transform with respect to [3 gives 

J dT e~°"/2(T, T) = exp (^^^^T - l^jTv/^M^) , 

(154) 

where putting a = Q confirms the normaUzation condition 
Jo°l2{T,T)dT = 1. Now by inverting the other Laplace 
transform with respect to a one gets the distribution 

/^(r, T) = ^£_e-M^(--i^)V2r^ (155) 



where t = t — T. 



X. INVERSE OCCUPATION TIME WITH DISORDER 

(<7>0) 

In the presence of disorder, i.e., a > in Eq. (|6j, taking the 
disorder average in Eq. ( I139l l gives 



/■oo /"OO 

/ dTe-P^ / dte-"*/occ(t|T) 



a+ p 



(156) 



where £2(a,p) is obtained by taking the disorder average of 
Eq. ( I108> . using the distributions of — z+(0) and z_ (0) given 
by Eqs. ( IB9> and ( IBl 1> respectively with a+ = a + p and 

a_ = a. 

It is useful to consider a different form of the above equa- 
tion, which is obtained by taking the disorder average of 
Eq. ( fT40l l 



dz e 



dr ( 









1) 



= 4(/3,a-/3), 



(157) 

where by taking the limit [3^0, one obtains the large T 
behavior of Iocc{t\T). 

A. Flat potential {fj. — 0) - Sinai model 



We will now study the large T behavior of Iocc{t\T), for 
the Sinai potential (/i — 0), by analyzing Eq. ( I157t in the 
limit /3 0. 

It follows from Eq. (|C3 that 



TOi(/3, a - /3) 



where mi (a, p) is given by Eq. iC3i . and 



(158) 



= 2/^0 and = 2ifo f ^ j • (159) 



In the hinit /3 ^ 0, since £7+ — log /3, hence 

-j-r — wi(0,a) 

|_-2i<ro (V2a/cr) log/3J 

which suggest the following scaling form at large T, 

1 



(160) 



locciAT) 



logT 



(161) 



Therefore in the limit /3 — > 0, substituting the above scaling 
form in Eq. (^^} and using Eq. ( I160> one gets 



/•oo 

Jo 



mi(0, a) 



2Xo (V2^/cr 



(162) 



However, the above Laplace transform is the same one given 
by Eq. jCl li in appendixIC] where a is replaced by p. There- 
fore we can directly borrow the results obtained there. Using 
the results from Eq. ( IC19> gives 



/3(r) « as T -> 0, 



and results from Eq. ( IC26t gives 



1 



^3(t) ~ — as r ^ cx). 
^ ^ 2t 



(163) 



(164) 



B. Unstable potential (/i > 0) 

For /i > 0, Eq. il5H suggests that in the large T limit. 



Ioccit\T) will only depend on the difference {t — T), 



Ioccit\T)^h{t~T). 



(165) 



Therefore in the limit /3 — ^ 0, using the above form in 
Eq. (I157> one gets 



dre-^^hiT) =4(0, a), 



(166) 



However, the above Laplace transform is the same one given 
by Eq. {D3\ where a is replaced by p. Therefore borrowing 
the results from appendixiDlreadilv gives 



h{T) 



(167) 



for small t and 

h{T) ^ e-^^ (168) 

for large r, with the same constant h as in Eq. ( ID16> . 

The normalization condition l4{T)dT = £i(0,0)=l/2, 
indicates that the particle escapes to — oo with probability 1 /2. 
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FIG. 10: The scaling functions go{x) in Eq. <170> plotted by using 
Eq. (m\ . 



C. Stable potential (fi < 0) 



We are interested in the behavior of Iocc{t\T) in the scaling 
limit where t ^ oo and T oo, but the ratio x — t/T is 
kept fixed. Substituting T — z/a,t — xz/a and p = as in 
Eq. (I156> . we get 



dx 



-(s+x)z 


-I ( 


xz 




-'occ 1 




-) 




a \ 


a 


a / 


= ilia 


, as) = 


1 


- ^2 



1 - m3(a, s) 

1 + s 
(169) 



Eq. ( IE6t . The above equation suggest the scaling form 

1 



/occ(i|T) = -5o(t/r), 



(170) 



with the normalization 9o{x) dx = 1, which follows di- 
rectly from the normalization Iocc{t\T) dt = 1. By substi- 
tuting the above scaling form in Eq. M69\ in the limit a ^ 0, 
after simplification one gets 



/"OO 


'x-l 


ll 


X + s 



go{x)dx = TO3(0,s), 



(171) 



where 7713(0, s) is given by Eq. ( IEI6I 1. By making a change of 
variable y = 1 /x, Eq. (IE16t reads 



7713(0,5) 



1 



x-l 



{x-iy 



■dx. (172) 



Therefore comparing Eq. ( I171> and Eq, M12\ readily gives the 
inverted Beta law 



9o{x) 



1 



{x-D 



x > 1, 



(173) 



which is displayed in Fig. [TO] The scaling function (7o(a^) has 
a maximum at x = 2v/{v + 1) for v > 1. However, go{x) 



diverges near x = \ for v < 1. Note that for v = 1/2, 
Eq. jl73l l gives identical results to that of a pure Brownian 
motion (/7 = and a — 0), given by Eq. J143t . 



XI. CONCLUDING REMARKS 

In this paper we have considered the motion of a particle 
in a one dimensional random potential. We have presented a 
general formalism for computing statistical properties of func- 
tionals and the inverse functionals of this process. We have 
used a backward Fokker-Planck equation approach to calcu- 
late the pdf of these functionals for each realization of the 
quenched random potential. The most difficult part of the 
problem is to carry out the disorder average on these pdfs. 
Thus to demonstrate the formalism explicitly, we have chosen 
the external potential to be the combination of a deterministic 
part and a random part, U{x) — ~'^\x\ + ^/aB{x), where 
B{x) is the trajectory of a Brownian motion in space. The 
case /i = in the potential corresponds to the Sinai model. 
The deterministic part of the external potential is stable for 
/i < and unstable for /i > 0. The pdfs of the functional 
and the inverse functional vary from one realization of B{x) 
to another, and in this paper we have shown how to carry out 
the disorder average on them, for two particular functionals, 
namely, the local time and the occupation time, and their in- 
verse. Despite the simplicity of the model, we get very rich 
and interesting behaviors by tuning the parameter /i/cr, which 
we have summarized in Tableslllllllandlllll for /i = 0, /i > 
and /i < respectively. In many cases the disorder changes 
the behavior of the pdf drastically from the pure case (cr = 0). 

A very interesting phase transition in the ergodicity of the 
particle position occurs at a critical value of the parameter 
|/i|/(T = 1, when the deterministic part of the potential is sta- 
ble (/i < 0). For I /i I/cr < 1, when the particle gets trapped 
in the wells of the random potential, the deterministic force 
— sign(a;) is not strong enough to lift it from the well and 
push it towards the origin and hence there are small number of 
zero crossings. On the other hand, for |/i|/CT > 1, the strong 
deterministic force sends the particle frequently towards the 
origin, and hence the system becomes ergodic. This change 
in the ergodic properties shows up in the qualitative change 
in the curvatures of the disorder averaged pdfs when the pa- 
rameter I' — I /i 1/(7 passes through unity. While for v <\, the 
disorder averaged pdf of the occupation time Pocc{T\t) is con- 
cave upward with a minimum at T = t/2 and diverges at both 
ends T — and T — t; for > 1 it is concave downward, 
which goes to zero at the two ends T ~ and T — t, and 
has a maximum at T = t/2 (see Fig.|8|i. In the context of in- 
verse occupation time, while for < 1, the disorder averaged 
pdf /occ {t\T) diverges near its lower end t = T and decreases 
monotonically as t increases, for z/ > 1 it has a maximum at 
t = \2v I (y + 1)]T and goes to zero at both ends t = T and 
t ^ 00 (see Fig.llO>. Similarly, the disorder averaged pdf of 
the local time Pioc{T\t) diverges near the lower end T = 
and decreases monotonically as T increases for v <\. On the 
other hand for v > has a maximum and goes to zero at 
both ends T = and T ^ 00 (see Fig.©. 
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For the stable potential, another very interesting observa- 
tion is that at |/i|/cr = 1/2, in the limit T ^ oo and t ^ oo 
while keeping the ratio T/t fixed, the exact asymptotic disor- 
der averaged pdfs of the occupation time Pace {T\t) and inverse 
occupation time Iocc{t\T) become exactly identical to the re- 
spective pdfs Pocc{T\t) and Iocc{t\T) for the simple Brownian 
motion (/i = and a = 0). It looks as if at the particular 
value \fJ,\/a — 1/2, the effect of disorder is exactly canceled 
by the deterministic stable potential. However, similar con- 
clusion is not true in the context of the local time and inverse 
local time. Therefore, a physical understanding of what ex- 
actly happens at this particular value of the parameter will be 
extremely useful. 

There are several directions open for pursuing research far- 
ther in this area. In this paper we have considered only the 
average of the pdfs over disorder. However, in many cases, as 
we have seen in this paper, the disorder broadens the distri- 
butions considerably. For example, for the unstable potential 
(/i > 0), even though for each realization of random poten- 
tial the local time has a narrow exponential distribution, by 
taking the disorder average one gets a broad power law dis- 
tribution, which is the indication of large sample to sample 
fluctuations and lack of self-averaging. Therefore, in this sit- 
uations the knowledge about the disorder averaged pdf (first 
moment) is not enough, and one requires to compute the other 
higher moments (over disorder). Thus extending our formal- 
ism to compute the full distribution (over disorder) of pdf will 
be very useful. 

The random part of the potential we have considered in this 
paper is very particular, where the barrier heights grow as -y/x. 
However, in realistic systems the random potential remains of 
order one throughout the sample. Therefore, it will be very 
interesting to extend this formalism for more realistic random 
potentials. 

Recently several asymptotically exact long time results for 
other quantities in Sinai model were obtained by using a real 
space renormalization group method ll39il . Using that method, 
reproducing the exact results obtained in this paper remains 
as challenging open problem. Another interesting direction is 
to study the properties of functionals of a more general non- 
Markovian stochastic process in random media, and to extend 
our results to higher dimensions. 



XII. APPENDIXES 

APPENDIX A: PDF OF THE LOCAL TIME IN THE CASE 
OF THE STABLE POTENTIAL, ^ < AND a = IN EQ. IS) 

In this appendix we will derive the pdf of the local time 
Pioc{T\t), for the stable potential (/x < 0) in the absence of 
disorder (a = 0). In this case by solving Eq. ( I23> with the 
boundary conditions oo) — > and ?/_(a; — > — oo) — > 

we get 



y±{x) — y±{0) e:xjp ^ ( — v + + 2a 



(Al) 



where ly = Substituting these results in Eq. i26\ we get 
A (a) = —1/ + \/v^ + 2a. Therefore the Laplace transform 
G[a) in Eq. ( l28l becomes 



■+ \/v'^ + 2a 



exp 



-V + Vi^^ + 2a j T . 

(A2) 



Now making a shift a — [3 — /2, in Eq. (I28> yields 



uT 







VP + v/^/2 



, (A3) 



where the right hand side is the Laplace transform of 
e" *'^^PiociT\t). Inverting the Laplace transform with respect 
to f3 and after simplification gives the exact distribution of the 
local time for all T and t. 



V2 



PUm - 4^e-(^-'^*)^/2'-^e^"'^erfc i^Vt+^], 



TTt \2 V2i. 

(A4) 

where erfc(a;) is the complementary error function. Note that 
Eq. (IA4> reduces to Eq. OH for v = 0. 
For large t, since 



, V r- T 

erfc ( — 

,V2 V2i 



X exp 



1 


.V2 


T 

V2i 


i- 


V2 


T ' 



Eq. (IA4> simplifies to 



V2 



-(T-i/t)^/2t 



TTt 



(A5) 



(A6) 



Putting = in the above equation one still recovers the 
result given by Eq. ( I31> . For non-zero ly, near the mean (T) = 
i^t, the pdf of the local time reduces to a Gaussian one 



Ploc{T\t) 



1 



-(T-ptY/2t 



(A7) 



APPENDIX B: PDF OF THE SLOPE VARIABLES z± (0), 
THAT APPEAR IN THE DISORDER AVERAGE 
COMPUTATIONS 

Both in the contexts of local and occupation time we have 
a homogeneous differential equation of the type 



1 



yUx)+F{x)y'±{x)-a±y±{x)^(), (Bl) 



with the boundary conditions y+(a; ^ oo) ^ and y-{x 
—oo) 0, and the force 



F{x) = /xsign(x) + V^^{x), 



(B2) 



with {£,{x)) — and (^(a;)^(x')) = 6{x — x'). For each real- 
ization of {C(a;)} in the force F{x), the solution of Eq. (IB1> 
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TABLE I: Flat potential. — Disorder averaged pdf 's of the local, inverse local, occupation and inverse occupation times of a particle starting 
at the origin, diffusing in the Sinai potential U{x) — ^/gB{x), where B[x) represents the trajectory of a Brownian motion in space with the 
initial condition -B(O) = 0. 



PURE CASE (o- = 0) 


DISORDERED CASE (a > 0) 




T/t fixed t log t 

y 




^ ' ' t/Tflxed Tlog^ y y V / J 

2 

gi{x) = - e^^^'^'^Ko (l/crx) 
a; 




Pocc(T|t) = i?L(rlf) + - Tjt) 

RL{T\t) ^ ^ i?(r) 

R{T) ^ for small T 
VTrT 

R{T)^^, for larger 




/3(r) ~ —z^, for small r 
/3(-r) ~ -!-, for larger 



is different, and for the disorder averaged computations per- 
formed in this paper we finally require the distributions of 
the stochastic variables y^. (0) /y± (0) and in this appendix our 
goal is to find these distributions. 
By defining the variables 



(B3) 



we find from Eq. JB1> that zj^{x\ satisfy the stochastic Ric- 
cati equation 



2a^ 



(B4) 



However, now the boundary conditions for z± in Eq. (I23> 
are not specified. Therefore for each realizations of {^(x)}, 
the solutions of z± (x) involve one unknown each that can not 
be eliminated due to the lack of the boundary conditions. In 
other words, to find the distributions of ^±(2;), we need the 
respective distributions at some initial points, unfortunately 
which are not specified. 

It turns out, however, that this difficulty can be bypassed 
by a method 11371 BSl l4lll which lets us to compute the dis- 
tributions of z+ (0) and Z- (0) without having the knowledge 
of the boundary conditions on z_|_(oo) and 00). We will 
present the method below for the present context. 



First we consider Eq. {B4\ for x > 0, i.e., 

= -zl{x) - 2[/i + ^/^£_{x)]z+{x) + 2a+. (B5) 

Note that z+(x) = y'_^{x)/y+{x) is negative. We make 
a change of variable x = — t and substitute z^{—t) — 
— exp[0(T)] in Eq. iB5\ to find that the new variable (/)(t) 
satisfies a much simpler stochastic differential equation 



rfr 



= 6(0) + 2V^e(T) 



(B6) 



where ^(r) = ^(-t) and thus (^(r)) = and (C(r)|(r')) 
6{t — t'). The source term 6(0) is given by 



2a. 



2fi. 



(B7) 



Now we can interpret Eq. (IB6> as a simple Langevin equa- 
tion describing the evolution of a Brownian particle starting 
at time r — —00, in a classical stable potential U(.i{4>) = 

- J^biifi) dip ^ e'l' + 2a+e-'> - 2/x0 - (2a+ + 1). Even 
though we do not know the starting position of the particle 
0(— 00), it is completely irrelevant. No matter what the ini- 
tial position is, eventually after a long time, i.e., when t is far 
away from —00, the system will reach equilibrium and hence 
the stationary probability distribution of is simply given by 
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TABLE II: Unstable potential. — Disorder averaged pdf of the local, inverse local, occupation and inverse occupation times of a particle 
starting at the origin, diffusing in the unstable random potential U{x) = —j-L\x\ + ^/aB[x), where fj, > Q and B[x) represents the trajectory 
of a Brownian motion in space with the initial condition -B(O) =0. We denote = ji/o. 



PURE CASE (cr = 0) 



DISORDERED CASE (a > 0) 



PiociT), Pioc(T')=2Ate 



-2iiT 



floc(T|t) 



floc(T), Pioe(T) =2^(l + ar) 



-(2i^+l) 



T 



V2^ 



exp 



2t 



+ (1 - exp[-2/ir]) 5{t - oo) 



52(2;) 



20F 



ir2(i/) 



+ {I - [1 + gT]-^") &{t - 00) 
^^^C/(1/2,1 + .,2M) 



Pocc(T|t) = RL{T\t) + RL(t- T\t) 
RL{T\t) Rl{T) 

exp 
1 



Pocc{T\t) = RL{T\t) + RL{t~T\t) 



T 



Rl{T) 



^expf^rUrfcf^^ 



/ttT 



2 'J \V2 
for small T 



RL{T\t) Rl{T) 

Rl{T) ^ , for small T 



/ttT 



for large T 



b is given by the zero of Ku{^/2p/o-) closest to origin in the left part 
of the complex-p plane. 



/occ(t|T) Ii(t -T) + h{t ~ 00) 



7i(r) =^i^/2 exp 
1 



/occ(t|r) J4(t - T) + - (X)) 

74 (r) ~ ^^^=, for small r 



AW 



/i\/2 



7rr 



7rr 

^/2 e-'-'"/^ 



for small r 



for large r 



74(1") ~ e , for large r 
fe is the same constant as above. 



the Gibbs measure 



where 



Pst(</>) = A exp 



= A exp 



(B8) 

where A is a normalization constant such that 
^st(0) ~ 1- Now changing back to the origi- 
nal variable zj^(x) we obtain the distribution of 2:+(0) as 



P+ (-z+(0) = w) 



-wf"^"-^ exp 



1 

2^ 



2a. 



(B9) 



w^/'^-i exp 



1 

2^ 





2a+r 




w J 



2(2a+)^/2^i^^ 



n/cr 



dw 



(BIO) 



Similarly for a; < 0, by putting F{x) — —fi + \/^(,{x) in 
Eq. (IB4> and substituting z_(a;) = exp[(/)(x)] one finds that 
(j){x) satisfies the same differential equation in x as Eq. (IB6> 
with (,{x) — —^{x) and a_|_ is replaced with a_. Therefore 
(f){x) has the same stationary distribution as Eq. ( IB8> and con- 
sequently the distribution of z_ (0) is same as that of — z+(0). 
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TABLE ni: Stable potential — Disorder averaged pdf's of the local, inverse local, occupation and inverse occupation times of a particle 
starting at the origin, diffusing in the stable random potential U{x) = — + ^/gB{x), where ^ < and B{x) represents the trajectory of 
a Brownian motion in space with the initial condition B{Q) = 0. We denote v = l/il/cr. 



PURE CASE (a = 0) 


DISORDERED CASE (a > 0) 


P\UT\t) ~ exp 
$(r-) = i (r 




J ' W/t fixed 
near r — /i 


/2(j/)= (^)"'"'e-'^/''t/(l/2,l + ^,2y/a) 

[(7i^(l^)J \(7/ 






53(2;)= .pjfy 7 ^-—.U{l/2,l + u,2/ax) 




2t \ 


Pocc(T\t) ~ exp 
$(r-) = 2^^ 


(- 


:?) 


1 ft ^ oo,T ^ oo 
J ' iT/t fixed 

1 

, near r = — 
2 


Pocc(T|t) jMT/t) 

T/t fixed t 

/o(y) = ;H7^[y(l-?/)]"'"'> 0<y<l 


Iocc{t\T)^l2{t-T,T)e{t 


-T) 

[ 2r J 


B(v,v) x''" 



namely 



p- (z_(0) = w) 



with 



0_ 



exp 



1 

'2^ 



2a_ 



(BID 



= 



1 

2^ 





2a_ r 




w J 



(B12) 



Note that the distributions of — z+(0) and z_(0) given 
by Eqs. ( IB9> and JBIU have maxima at (/i — <t) + 

— fj)^ + 2a± respectively and in the limit a 
the distributions tend to delta functions around their 
maxima. Therefore in the limit a Q one re- 

covers the pure case results by using the distribu- 



tions P+ (z+(0)) = 5 [z+{Q) + [m + _^ 2a+] j and 
P- (z_ (0)) = 5 (z_(0) -[^Ji+ v//^2 + 2a_]) . 



APPENDIX C: LEFT HALF OF THE DISORDER 
AVERAGED PDF OF THE OCCUPATION TIME FOR SINAI 
POTENTIAL = AND ct > 0) 



By taking the disorder average of Eq. \\ \2\ one gets 

/ dte"^' dTe-P^ RL{T\t) ^ -ei{a,p). (CI) 
Jo Jo Oi 

Using the distributions of — z+(0) and z-{Q) from Eqs. (IB9> 
and (IBIU respectively with a+ = a + p and a_ = a, from 
Eq. ( I107t one gets 



£i(a,p) = 



™i(q;,p) 



(C2) 



where, 

mi(Q;,p) 

X 

and 



dull 
wi 

dW2 

W1+W2 



exp 
exp 



1 

'2^ 



Wl 



1 



2a 



2cr V W2 



12^ 



(C3) 



(C4) 
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Before we proceed further, let us take a detour to check the 
normalization condition of i?x,(T|t). By putting p = in the 
above equations we get 



and 



■mi{a, 0) 



n^ = n = 2Ko 



°° r°° dwi dw2 

Jo Wi W2 

X exp 
X exp 



/2a 



W2 



(C5) 



Wi + W2 



1 




2a 




wi H 




2^ 






1 




2a 




W2 H 




2^ 




W2 



(C6) 



Note that the above integral must remain invariant under the 
transformation wi W2 of the dummy variables. Therefore 
we get 



2mi(a,0) 



' dw r 1 




2a' 




/ — exp<^- — 


w - 




}] 


.Jo ^ I 2(7 




w 





(C7) 

Therefore we have £i{a,0) — 1/2, and inverting the Laplace 
transform in Eq. (ICU with respect to a for p = gives the 
normalization condition dT = 1/2. 



Now we analyze the large t behavior of i?o(r|t). By mak- 
ing a change of variable z = at, it follows from Eq. (ICU 



/ dze"^ / dTe-P'^RoiT,z/a) ^ ei{a,p). (C8) 
Jo Jo 

In the a ^ limit, 51+ 2Kq [y/2p/a) and il_ - log a. 
Therefore from Eq. ( IC2l l 



mi(0,p) 



-2ifo(V2^/^) log a] 



(C9) 



which suggest the following form for RL{T\t) at large t 

1 



RLiT\t) 



\ogt 



RiT), 



(CIO) 



where R{T) is independent of 

Now using Eqs. ( IC9l l and dClOt . in the limit a ^ 0, 
Eq. iCSl gives 



/ dTe-P'^R{T)^ 
Jq 



TOl(0,p) 



2i^o (V2^/fT) 
where mi (0, p) is obtained from Eq. (IC3t . 



(Cll) 









/ exp 






/o ■"'1 




WlJ 



X / c?it;2 • 

Iq Wl + W2 



(C12) 



By making change of variables wi = 2ax and W2 = 2ay in 
the integrals in the above equation one gets 



mi(0,p) 



dx 




/ — exp 


[ (^+2.2 J] 


/o 2^ 





/o 2/ + a; 



(C13) 



where now the integral over y can be expressed in terms of the 
incomplete Gamma function i44ll 



dy — ■ — 

y + x 



e^r(0,x) 



(C14) 



Therefore, after straightforward simplification, Eq. ( IC13> be- 
comes 



mi(0,p) 



P°° dr 

/ —e-'=ri0,p/2a^x). (C15) 
Jo ^ 



Now we will analyze the limiting behavior of R{T) for 
small and large T by taking the limit of large and small p 
respectively. 

Since for large p 



T{0,p/2a^x)^^^. 

P 

from Eq. ( IC15> , one gets 

2o-2 1-°° 

mi(0,p) w / dx exp 

P Jo 

2V2a 



-p/2a^x 



(C16) 



(^+2^) 



VP 



(C17) 



Since the asymptotic behavior of ^ (x) is independent of 
substituting ?7ii (0, p) from above in Eq. dClU gives 



/"OO 

/ dTe-P'^R{T) 
Jo 



A/2cr 



(C18) 



for small p, and by inverting the Laplace transform with re- 
spect to p one obtains 



R{T) 



V2a 

7^' 



as T 0. 



(C19) 



To obtain the large T behavior, we first consider the follow- 
ing integral 

r°° dr 

V{z)^ — e-=^"r(0,p/2(j2a;), (C20) 
Jo X 

where X'(l) = mi(0,p), follows from Eq. ( IC15> . Now by 
differentiating ^{z) with respect to z, one can express it in 
terms of the modified Bessel function as |44] 

V'{z) = - dxe-'='r{0,p/2<j^x) ^ ^^Ko (^^^) • 

(C21) 
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Now by integrating back again with respect to z, we obtain 

mi(0,p) as 

mi(0,p)=I?(l)=2^°°^ifo(^) 



f°° dx 
4 1 —Ko{x), (C22) 

/2p/(T ^ 



where we have made the change of variable 2pz/a'^ — x^ . 
The p ^ Hmit can be obtained from the hmiting behavior 
of the integral 1 45 1 



KQ{x)dx 1 , .2 „ 

;^(logy) asy->0, (C23) 

X I 



which gives 



1 2 

"^i(0,p) ~ - (logp) , asp-^0. 



(C24) 



as p ^ 0. Thus, inverting the Laplace transform with respect 
to p one obtains 



KT)-^. asT^oo. 



(C26) 



APPENDIX D: LEFT HALF OF THE DISORDER 
AVERAGED PDF OF THE OCCUPATION TIME FOR /i > 



By taking the disorder average of Eq. W 12t one gets 

/ dte-"' dTe-P'^ RL{T\t) ^ -ei{a,p). (Dl) 
Jo Jo a 

As in the pure case (a = 0), one also expects the t ^ cx 
behavior of the disorder averaged distribution to tend to a i 
independent form 



t — >oo 

Therefore Eq. (IDU becomes 



lim RLiT\t) = Rl{T). 



(D2) 



/ dTe'P^RjT)^h{Q~pj. (D3) 
Jo 

Using the distributions of — z+(0) and z-(0) from Eqs. (IB9> 
and (IBl 1> respectively with a+ = p and a_ = 0, from 
Eq. (I107> one gets 



1 f°° 



exp 



— ^ — 



2(7 



X / dW2 — . 

/o W1+W2 



(D4) 



(D5) 



with 

n+ = 2(2p)^/2-x^/, 

= (2a)'^/'"r(/i/(T) 
Now the integral over W2 in Eq. ( ID4> can be expressed as |4; 



dwo ^ 



Wi + W2 

w 



./V"./2-rf^ + i)rf-^,^V (D6) 

Vcr / \ a 2a ) 



where r(a, x) is the incomplete Gamma function. Therefore 
Eq. (ID4> becomes 



^1(0, p) 



[im2{p) 



(2a)A'/-+i (2p)A'/2-/f^/,(y2^/a) 



(D7) 



Since {\/2p/a) ~ — i logp, as p ^ 0, Eq. jClU gives where 

1 

dre-f^i?(T)^--logp, (C25) 



The small and large T behavior of Rl{T) can be found by 
analyzing Eq. \D1\ in the limiting cases p ^ 00 and p ^ Q 
respectively. 

Making a change of variable wi ~ p/ax in Eq. (ID8> . and 
then taking the p —> 00 in the incomplete Gamma function, 

gives 



m2( 



(p) « (2a)2(2p)^/^-iy da; x-^/'^ exp [- (a^ + ^) 

(DIO) 

where the integral above on the right hand side can further be 
expressed in terms of the modified Bessel function as 1I44I1 



dxx-^'/'' 
= 2{2aY''''\2pr^^'^-+''^K^,,_^ 









exp 




2cr^xJ\ 



(DID 



Since the large x behavior or K^{x) is independent of ly, 
Eq. (ID7> simplifies to 



Zin),p) « — — as p ^ 00. 



(D12) 



Therefore by inverting the Laplace transform in Eq. {D3t with 
respect to p, one gets 



Rl (T) w , for small T. 



/ttT 



(D13) 
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Now we will analyze the the large T behavior, by taking the 
limit p — > in Eq. jD7> . It is straightforward to obtain from 
Eq. (inD that 



m2(0) = 



(D14) 



which gives 



(D15) 



for small p. However if one takes the Umit p now in 
K^/^{^/2p/a) in the above expression, it only gives the nor- 
malization condition Rl{T) dT — 1/2 and does not pro- 
vide any information about the large T behavior of Rl [T). 
We make the following ansatz 



Rl{T) 



for large T. Then the Laplace transform 

r°° 1 

/ dT e^P'^ Rl{T) K — 
■/o P + 



(D16) 



(D17) 



for small p. Therefore substituting Eqs. ( ID15t and ( ID17> 
m Eq. (|d3, one can conclude that b is given by the zero 
of K^/„{^j2p/a) closest to origin in the left part of the 
complex-p plane. 



APPENDIX E: DISORDER AVERAGED PDF OF THE 
OCCUPATION TIME FOR ^ < 

Taking the disorder average of Eq. ( fT06l 2ives 

/ die""* / dT e~P'^ P„cc{T\t) 
Jq Jo 



a+p a{a + p) 



£i{a,p), (El) 



where we have substituted £2{a,p) = 1 — £i{a,p). 



We are interested in finding the behavior of Pocc{T\t), in 
the scaling limit t oo,T oo, but keeping T/t = y 
fixed, which corresponds to the limit of conjugate variables: 
a — > 0,p ^ 0, keeping p/a — s fixed. Substituting z = at, 
T = yz/a and p = as in Eq, iEl\ . we get 



dy 



dze-^^+'y^' 



Ip (Ml 

^ occ 1 

a \ a 



_ 1 + 57713(0;, s) 

(1 + s) ^ 
(E2) 

where 7773 (a, s) = ii{a,as). Equation (IE2> suggest the form. 



Pocc(T|t) = - UT/t), 



(E3) 



in the scaling limit t ^ 00. T — > 00, while their ratio T/t 
is kept fixed. In the limit a ^ 0, by substituting the above 
scaling form in Eq. JE2> . it is straightforward to obtain 



' foiy) 1 + 57773(0,5) 

dy ■ - 







1 + sy 



+ 



(E4) 



where putting s — 0, gives the normahzation condition 
Jq fo{y)dy — 1. Using this normalization, the above equa- 
tion can be simplified to the following elegant form 



' i-y 

l + sy 



fo{y) dy = 7773(0,5). 



(E5) 



Now by using the distributions of -2;±(0) from Eqs. ( IB9l l 
and (IBl It with a- — a and a+ = a{l + s), from Eq. J107t 
we get 



7773(0, 5) = £i{a, as) 
1 



X w. 



dwi 



exp 



dwo 



W2 







X w. 



exp 
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'2^ 
1 



Wi + W2 

2a(l + s) 



Wl 



Wi 



2a 



„ , W2 H 

2ct I W2 



where 



= 2(2a)-''/2(l + s)-^/'K, ( + 



(E6) 

(E7) 
(E8) 



with 1/ = l/^l/f . Note that we simply can not take the limit 
a — > in the integrals in Eq. iE6\ . as it diverges in that limit. 
However, it is possible to extract the divergent contribution 
outside the integrals which finally cancels exactly with the di- 
vergence of n±. This is done by making the change of vari- 
ables 



a(l + 5) 
in the integral to get 



7773(0, 5) 



X, and 



(1 + .'^)- 



a 



aw2 



y, 



(E9) 



dx / dy 
Jo 



a; + (1 + s)y 



X exp 
X exp 



a(l + 5) 
2a'^x 
a 



y 



2a^y 



(ElO) 



Now the limit o ^ can be taken in the above equation, 
as, in this limit f7+ -> cr''a~''(l + s)^T{iy) and 



a'^a Therefore from Eq. ( IE1Q> we get 

"^3(0,s) = / dyy^-^e^y —————dx. 

W Jo Jo x + {l + s)y 

(Ell) 

Now the integration over x can be expressed in terms of the 
incomplete Gamma function 14411 

-ty-p l-oo -X 

r(p,A) = — ——dx, [Rep<l,A>0], 

(E12) 

which gives 

m3(0, s) = \J > / y2.-ig.,r (1 + s)y) dy. 

(El 3) 

The right hand side, however, is one of the inte- 
gral representation of the Gauss's hypergeometric function 
F{a, P; 7; z) 13, which gives 

m3{0,s) = ^Fil,iy;2iy+l;~s). (E14) 
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